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TECHNICAL MEMORANDUM 


SOLAR ARRAY FLIGHT EXPERIMENT 
BACKGROUND AND PURPOSE 


The purpose of tliis study was to predict the longitudinal frequencies of the solar blanket for all 
states of deployment. Furthermore, a study was done to detect a possible crosscoupling between bending 
and longitudinal frequencies, Finally, the study attempted to predict potential problems caused by the 
station keeping maneuver. For better understanding of this problem, Figures 1 and 2 were included to 
show the construction of the SAFE hardware. 


DESCRIPTION OF THE SAFE STRUCTURE 


The SAFE structure consists of a deployable mast and 84 solar panels that are liinged sequentially 
on torsional springs, forming a deployable blanket. One end of the blanket is attached to the header of 
the mast; the other end is connected to the base of the mast. 

In this study, an assumption was made that all boundaries around the blanket were rigid. At all 
times, the solar panels were guided by wires to deploy simultaneously with the mast. 

The motion of the blanket is restricted to the.xy-plane, consequently, a sail-type bulging is 
excluded (Fig, 3), 

The first section consists of 58 panels or an equivalent of a system with 28 degrees-of-freedom. 
Deployed, it represents 70 percent of the total length. 

For computation, the above system is restricted to six panels only (Fig. 4). 

During the deployment each panel rotates and translates simultaneously. For computation, the 
rotational points chosen are as the center of mass of each panel. The motion of each panel is fully 

described by the angles 0j'. Consequently, the angles become the “generalized coordinates.” The 
geometry of the system requires that the angles are; 



( 1 ) 


Only two points (Pj and P 2 ) on the six-panel blanket of Figure 4 are free to move anywhere 
within the boundaries. This narrows down the degrees-of-freedom to 





( 2 ) 



N = degrccs-of-freedom 
n/2 = panel pair . 


THE MATH MODEL 

All variables are expressed in natural coordinates Q(, and the constants are: 

n = 84 total number of panels 

1 = 15 in. width of the panel 

k = 0.795 in. Ib/rad torsional stiffness of the liinge 

m = 0.0106 lb sec^/in. mass of a single panel 

I = m/3 • (1/2)^ = 0.1992 in. lb sec^ inertia mass of a single panel 

= k/I = 3.992 rad/sec^ = a chosen abbreviation, called “the reference frequency" 
cOq a 2 rad/sec 

kinetic energy 


T = 





NOTE: [see equation (6)]. 

Potential energy without gravitational contribution. 


V = |(6i')2 + |2(6i_,' 


i=l 


+ 0i' )^ + j (^n' 


Lagrange’s equations: 



(3) 


(4) 

h 


( 5 ) 


2 


The equations of motion are obtained by substituting equations (3) and (4) into equation (5), 

To keep the derivation simpie, a math model with 2 degrees-of-freedom is used. Later, tlie 2 degrees-of- 
freedom matrices arc extrapolated to represent 28 degrccs-of-freedom. This tcclinique is applicable only 
if the matrices arc symmetric. Because of Lagrange’s equation property, the mass and stiffness matrices 
arc symmetrical. 

Tiic derivation is cumbersome but straight-forward. Using Figure 4, tlie center of mass of each 
panel can be dcscr‘’'ed with trigonometric functions using a common origin. By differentiating the 
displacements of each panel with respect to time, the linear velocity (vp is established, and the rotational 

velocity needed in equation (3) is a time derivative of the angles between the panels 6-y In order to keep 

tlie differential equations solvable, the assumption was made that the angles between the panels consist 
of a constant part and a time-varying part. Tiiis assumption implies that tlie panels oscillate around the 
equilibrium point 0^ 

The above assumption is defined as: 


Oil(t) = Oo + 0i(t) 


( 6 ) 


where 


= total angle between the panels 
0Q = angle of deployment or the equilibrium position 
Oj = oscillatory angle 


Using Figure 4, the displacement of the center of mass for the first panel is; 


xi = sindj^ =^1") sin (00 + di) 


(7) 


where the time derivative is: 


■ 6 )‘ 


Xl = I r jdl cos (0Q + dj) 


for 


lim 

d J-+0 


Xi = (l) dj cos do corollary ^ sin d^ 


3 


Consequently, the total linear velocity of the first panel’s CM is: 




^ cos^ + sin^ 0 q)1 


For simplicity, s = sin^ Oq, c = cos^ <?q. The rest of the velocities can be easily seen from Figure 4. 





r) 1(3^ 1 ) 


r) f(40r 


i) f(40i- 


3 \2 . . 

-j [(401- 


3\2 . 

r f(40i- 


2 c + s] 


[(401 + 30*3)2 c + 03^ s] 

[(401 + 403 + 05)2 c + (05^ s] 
[(401 + 403 + 305)2 c + 05^ s] 


1st panel pair 


2nd panel pair 


3rd panel pair 


( 8 ) 


NOTE: Because of equation (1), only odd or only even subscripts are needed. If one substitutes equa- 
tion (8) into equation (3), the kinetic energy is: 


■?(:)' 


0l2c 

+ «l2 

+(30j)2 c 

+ 6,2 

+(401 +03)2 c 

+ 632 

+(401 + 303)2 c 

+ 632 

+(40 1 + 40 3 + 0 5)2 c 

+ 652 

+ (401 + 403 + 305)2 c 



+ 7 (01^ + 032 + 052)*2 


( 9 ) 


4 


Because of the rigid boundaries, the end requirements arc: 


n n n 

V Ojl =53 (^0 + <?i) = n % = constant 

i=I i=I i=l 


or 

n 

53 “ constant - hOq = 0 

i=l 

in this case 


<?i + O 2 ^6 “ ® 


since 


or 


Oj = (?2 » 


0 j+ 03 +Og — 0 O5=-Oj-03 , 

The mass moment of each panel rotated around CMp^^j^gL inertia of each panel is. 

"T-ft)’ ■ 

Using the above restrictions, the following substitutions were made: 


T = |l 


01^c 


+ 0J^ s 


+(30i)2c 


+ 0 s 


+( 40 1 + 03)^ c 


+ 0323 


+( 40 ’j + 303)^ c 


+ 032 s 


+(40 j + 403 “ 0 j 

- «3)2 0 

+ ( 0 i 3 + 03 )' 

+(40 j + 403 - 30 

- 393)2 c 

! + (0 1 + 

03)2 


+ I[ 0 j 2 + 03 ’^ + ( 0 j + 03 )^] 


( 10 ) 


( 11 ) 


( 12 ) 



; Using the cliain rule on equation (12) and differentiating the kinetic energy partially in respect to 

Oj and later witli respect to time, the equation of motion becomes: 


^ (Uj) “ ^ {3(480i + 24 O 3 ) cos^ Oq + \6d I + 

^ (If ) ° I {3(24ff| + I 0 O 3 ) cos^ + 85, + ISSjI 


Using Figure 4 and equation (4) the potential energy is; 


V = |(«r)2 + |(0l' +02')^ + |«'2' + +»4’)^ + |t»4' +«5') + |ces' 




because of equations (1), (6) and (10), V reduces to: 




V = 2 + ^1) + 2 + 2 + <^1 + ^3) + i (2 £>o + 203) + J ( 20„ - 0J) 

+ |(2«o-2»l -2«3)2 + |«I„-()| - 43 ) 


according to equation (5) the partial potential energy is: 


^ = k(0+ 1201 +603) 


3V 

900 


= k (0Q + 601 + 1003) 


If one collects all the terms, the final equations of motion are: 
r|3(480i + 2403) cos^ 0Q + I60’i + 80^^ + k(0 + 120J + 603) = 0 

I |3(240i + 160*3) "** ^®1 ^^^3} "** ^(®o 10^3) 0 


(13) 


6 



Abbreviatcti for convenience, in short: 



( 14 ) 


If equation (14) is substituted In equation (13), the equations of motion are; 



+ 3 cos'^ Q, 


48 24l j 1 2 r ‘2 6 1 ( 

24 16 1 ^3 j 6 I0_ i 



(15) 


Using symmetry, the equation of motion for a system with 5 degrees of freedom is; 
(MJ (0|)+VtKl {0i) = -w„2 


The equation is explained as follows: 



[M] = mass matrix consists of two parts, static and ..jrnamic. 

[K] = stiffness matrix 
* = multiplication sign. 

7 



DISCUSSION 


Equation of motion (15) was trefi?»^d as nn eigenvalue problem. A small desk computer was used 
to compute tlie roots of the problem since equation (15) can be rewritten as: 


IM] = Xi [K] 


or 


In Figure 5 the natural frequency was drawn as a function of tlie deployment angle 

As expected, the natural frequency of the blanket increases shortly before full deployment because the struc- 
ture becomes stiffen At full deflection, all the blanket modes (longitudinal) coincide at f ss 0.3 Hz. 


Using the finite element method, the dynamics of the wliole system were simulated. For 70 
percent deployment, the natural frequency of it is 0.06 Hz (Fig. 5, point R 2 ). To answer the question 

of a potential problem, due to a mutual excitation between blanket and mast, the frequency response 
curve ( 2 ) was drawn. Points Rj and R 2 were given by simulation and R 3 was implicitly computed using 
the beam equation. 

The curves (1) (Fig. 5) intersect with curve (2) at Sj, S 2 , S 3 , ... S 2 g. Those are points of mutual 

coupling that if given sufficient time, can build up to a problem - causing resonance. To investigate the 
time the mast is exposed to the sympathetic oscillation of the blanket, the dwell-time concept was 
introduced. 


DWELL FACTOR Tj„ 

T^^ is a factor that defines the percent of the mast period, the mast and blanket are essentially 

at coincident frequencies. Tliis factor is defined as the ratio, in percent, of the time a changing blanket 
natural frequency lies between the 1/2 power points of the mast first mode response, and is shown by 
the formula on Table 1. As shown in this table, the lower blanket modes (i = I, 2, 3) coincide with the 
mast first natural frequency at Sj, S 2 , and S 3 with frequencies fj = 0.070, 0.082 and 0.092, Since the 

blanket natural frequency is only coincident with the mast natural frequency for less than 1/4 of a 
mast period, the response does not have sufficient time to grow to a significant level. 


CONCLUSION 


An obvious problem is not evident. A possible further study could determine the dynamic 
behavior under the station-keeping conditions. 

S 


TABLE 1. THE DWELL FACTOR FOR THE LOWEST THiREE MODES 



Sj = rate of change of blanket mode frequency at tlic intersection with the 
mast first bending mode (Fig. 5). 

fj = frequency at intersection 

i “ t modal number 

j 0.0 1 

T= I/fj 

AT = 2f f;2/S| 

Tj,v - AT/T * 100% 










Figure 3. The overall view of the SAFE model. 



Figure 4. A closeup section of the blanket. 
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